Multimer formation in ID two-component gases and trimer phase in the asymmetric attractive 

Hubbard model 
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We consider two-component one-dimensional quantum gases at special imbalanced commensurabilities 
which lead to the formation of multimer (multi-particle bound-states) as the dominant order parameter. Lut- 
tinger liquid theory supports a mode-locking mechanism in which mass (or velocity) asymmetry is identified as 
the key ingredient to stabilize such states. While the scenario is valid both in the continuum and on a lattice, the 
effects of umklapp terms relevant for densities commensurate with the lattice spacing are also mentioned. These 
ideas are illustrated and confronted with the physics of the asymmetric (mass-imbalanced) fermionic Hubbard 
model with attractive interactions and densities such that a trimer phase can be stabilized. Phase diagrams are 
computed using density-matrix renormalization group techniques, showing the important role of the total den- 
sity in achieving the novel phase. The effective physics of the trimer gas is as well studied. Lastly, the effect of 
a parabolic confinement and the emergence of a crystal phase of trimers are briefly addressed. This model has 
connections with the physics of imbalanced two-component fermionic gases and Bose-Fermi mixtures as the 
latter gives a good phenomenological description of the numerics in the strong-coupling regime. 
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I. INTRODUCTION 

The notions of quantum liquids and their instabilities are 
paradigmatic for condensed matter physics [1]. For multi- 
component fluids, an important set of instabilities is associ- 
ated with interactions between components. A classic exam- 
ple is the Cooper instability of a spin- 1/2 Fermi liquid : even 
an infinitesimal attractive coupling between fermions of oppo- 
site spins drives a phase transition into the Bardeen-Cooper- 
Schrieffer superconductor li2J. A one-dimensional (ID) coun- 
terpart of the Fermi liquid, the spinful Luttinger liquid, has 
a similar instability, where an attractive inter-spin coupling 
opens a gap in the spin channel IUHl- 

Traditionally, the bulk of the discussion on two-species liq- 
uids assumed the SU(2) spin symmetry. The recent years 
have witnessed a growing availability of expeiimental studies 
of mixtures of unlike particles. This includes loading ultra- 
cold atoms to spin-dependent optical lattices |5], and trapping 
atoms of different masses |6] or even different statistics 17|. 
While most of experimental progress so far is in the domain 
of ultracold atoms, we stress that the relevance of such asym- 
metric mixtures is not confined to the realm of cold gases : 
dealing with more traditional solid state systems, one faces an 
asymmetric mixture situation as soon as the Fermi level spans 
several bands (which a priori need not be equivalent). This 
setup is typical for such diverse materials as semi-metallic 
compounds, mixed-valence materials, organic superconduc- 
tors lUt], small radius nanotubes and even graphene-based 
heterostructures [flO^. 

A generic question immediately arises : given a two- 
component mixture, what is the role of (the lack of) SU(2) 



symmetry? Or, more precisely, does the symmetry between 
components limit the set of instabilities of a liquid? Clearly, 
the answer might depend on the universality class of the liq- 
uid, and on the particular way the symmetry is broken. The 
simplest albeit non-trivial way of breaking the symmetry is to 
assume species-dependent masses of the particles. Even if we 
consider the few-body problem, this is known to bring new 
physics like the Efimov phenomenon : while for an equal- 
mass Fermi liquid the only allowed bound state is a Cooper 
pair, three -body bound states (trimers) appear once the mass 
ratio exceeds a certain threshold ifTH . The atom-dimer scatter- 
ing is strongly affected by the mass asymmetry fT?) and the 
ultimate fate of a Fermi liquid in presence of the Efimov effect 
is cuiTently an open question being actively investigated ifTsll . 
All these theoretical considerations are strongly motivated 
by cold-atom experiments which have recently achieved de- 
generacy of Fermi gases with different masses |j6|] and spin- 
imbalanced two-component fermionic gases [iTill . 

The physics of ID quantum many-body systems offers 
powerful methods f7], both analytical and numeiical, to have 
quantitative predictions on the fate of the Luttinger liquid in 
the presence of perturbations. The role of mass asymmetry 
for a two-component Luttinger liquid has been investigated in 
the renormalization group (RG) framework originally in the 
context of soUd state physics ifsi ITsi [T6ll . and recently revis- 
ited mostly in the context of cold atoms lfT7l - l25ll . and supple- 
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FIG. 1 . Qualitative formation of trimers in the weak-coupling picture 
of the bosonization (on the left) and in the strong-coupling picture for 
large interactions (on the right). 
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men ted by numerical investigations 11261 - 12911 . Overall the con- 
sensus was that the only new instability arising due to asym- 
metry is the collapse (demising) instability for large asymme- 
try and/or strong interspecies attraction (repulsion). Recently, 
a novel family of instabilities was predicted [30] to exist due 
to the interplay between polarization and asymmetry : These 
instabilities only take place for polarized mixtures of either 
statistics, and are characterized by the locking of the ratio of 
the densities to a rational value. Subsequent work in Ref. 3 1 
elucidated the relation of these instabilities and existence of 
few-body bound states. A qualitative picture of the mode- 
locking mechanism and the strong-coupling limit of the trimer 
formation is given in Fig. [T] The latter regime recalls an- 
other approach to multi-particle bound-states which is the use 
of many-colors (A^-component) fermions [32- 34] with which 
the physics of the trimers share qualitative features. 

This paper is divided in two main parts : the first one in- 
vestigates in detail the bosonization approach and the mode- 
locking mechanism mentioned above, while the second is ded- 
icated to the specific but important example of the ID asym- 
metric Hubbard model using the density-matrix renormaliza- 
tion group (DMRG) technique ||35ll . The predictions of the 
first part account for most of the numerical data, but a more 
phenomenological Bose-Fermi picture is proposed as a com- 
plementary analysis. Other important questions such as the ef- 
fect of a trapping potential or the emergence of crystal phases 
are eventually addressed. 



II. BOSONIZATION ANALYSIS 

In this section, we describe the salient features of the effec- 
tive bosonic field theory appropriate to a ID mixture of two 
distinct fermionic (or bosonic) atoms. The aim of this section 
is to give a bosonization interpretation for the formation of 
few-body bound states and their effective behavior through 
a mode-locking mechanism between the two species. Pre- 
dictions on the nature of the resulting phase are then made. 
The theory is a priori valid for models in the continuum or 
the continuous version of lattice models at generic (i.e. non- 
commensurate) densities. The effects of the presence of the 
lattice on certain commensurate densities will be briefly dis- 
cussed in section HTDI Notation conventions are standard and 
taken from Ref. |4j. 



A. Mode locking mechanism 

The two species are labeled by a pseudo-spin index a — 
t, 4- and their corresponding densities ricr such that n = + 
ni is the total density. Each species can be described by a 
scalar field (j)a- and its dual 9a- The creation operators can be 
expressed as a function of these fields, with, for fermions 



(1) 



and, for bosons. 



,i2p{k„x-(t>„) -i0 



(2) 



We have included all higher harmonics : as a consequence, the 
summation is over all integers fi^. The "Fermi momenta", 
ka — Trrig. are a priori not equal to each other, corresponding 
to a spin-imbalanced situation. The density operators read 



(3) 



The effective low-energy Hamiltonian can be written in terms 
of the fields (p^ and their canonically conjugate momentum 
Ho- = dxda/TT- In the case of absence of inter-species in- 
teractions the effective bosonic theory is given by Ho{<f>f) + 
'Hoih), where 



-Hoi 



2tt 



dx 



(4) 



where is the sound velocity and the so-called Lut- 
tinger parameter, which is equal to one in the free fermions 
or free hard-core bosons cases. Taking into account density - 
density interactions between species, of the generic form 
Jdxdx' U{x — x')h^{x)ni{x'), changes the effective theory 
and brings new kinds of terms : zero-momentum terms in 
the density representation Q couples the two spin species 
through a bilinear operator 



Hi 



9_ 
27r 



dx {dx(l)^){dx(f>i) 



(5) 



where g is a forward scattering constant, and higher harmonics 
terms involving multiples of the spatial frequencies : 

E ^pg / [27r(pn^ ± qni)x ~ 2{p(f)f ± gt/);)] (6) 



p,<J>0 



where G^^ are non-universal coupling constants. Clearly, if 
the generalized commensurability condition 



pn^ - qni ^ 



(7) 



is satisfied (with p and q coprime integers), and provided these 
terms are relevant, they will tend to lock the up and down 
fields together When the densities are fine-tuned to the def- 
inite commensurability (|7]i, then all other cosine operators in 
the sum are oscillating, in which case they don't contribute in 
the continuum limit (or they are less relevant for multiples of 
p and q). The remaining important operator in the sum (|6) is 
thus the sine-Gordon term 



Hi — G I dx cos 



with the combination 



72 



(8) 



(9) 
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For attractive interactions G < (we will argue below that 
this choice favors the relevance of the term), energy will be 
minimized when the field is pinned to (0a) = 0. Notice again 
that the above argument on the mode-locking mechanism does 
not rely on the presence of a lattice. Lastly, the cosine locks a 
combination of the bosonic modes but at a generic total den- 
sity n, there remains another bosonic mode leaving the full ex- 
citation spectrum gapless. We will see that the latter describes 
the effective behavior of the bound-states. In the following, 
we dub (pb this massless bosonic mode. 

We can draw a last remark on the operators (|6} : they have 
high scaling dimensions near the free fermion fixed point and 
are expected to be irrelevant apart from some special circum- 
stances which are the object of this work. In the fermionic 
language, they involve {p + g)-body interactions of the form: 

Eri^k(^o^ii(fco nvit(^;)^«t(fcf)+h.c. , 

{k} 1=1 j=l 

(10) 

where the summation over {k} runs over all combinations of 
2{p+q) — l momenta due to the total momentum conservation 
law : EUiih - K) + - ^f) = 0. Such interac- 

tions appear at high order in perturbation theory in a Hubbard 
model for example or after several steps of a RG treatment. 
For practical purposes it is simpler to work with the bosonic 
formulation given by Eq. and this is what we do from now 
on. 

In the following, we assume that the densities are commen- 
surate via the condition (|7]i, and analyze the simplified effec- 
tive theory written in terms of the t and | fields 



where the velocities v^r and Luttinger parameters are de- 
termined by the intra-species interactions. The quadratic part 
T~(-o{4>t) + '^oi4'i) + "Hi can be diagonalized by a Bogoh- 
ubov transformation |37] which could give a starting point for 
a perturbative RG calculation S [13 El IHQ . Due to the 
velocity asymmetry, additional couplings are generated and 
velocities are renormalized along the flow. The discussion 
of the nature of the gapped phases and their correlations re- 
mains unclear. In particular, diagonalizing the quadratic part 
T^oi't't) + ^o(04,) + ^1 of the Hamiltonian ( fTTl i does not 
give, apart from special choice of the parameters, the com- 
bination ^ that appears in the cosine term H2- In the next 
section, we take the following strategy : we look for the con- 
ditions under which the quadratic part and the cosine term are 
simultaneously diagonalizable. At the price of a restriction on 
the parameters range, the analysis can be done safely both for 
the criteria of relevance of the cosine and for the correlation 
functions in the single-mode phase. In spite of the limitation 
of the approach, we believe the scenario does occur without 
this restriction : as shown numerically in Sec. |lll] on a real- 
istic model, the single-mode multimer phase can span a wide 
region of the phase diagram. 

Lastly, we notice that, similarly to the phase separation 
criteria in two-component mixtures (when one of the mode 
velocities vanishes), the single-mode phase will undergo a 



phase separation instability when the gapless mode velocity 
Vb vanishes. We thus expect to find the single-mode phase 
surrounded with the two-mode phase and a demixed phase. 



B. Field transformation 

We have qualitatively discussed the fact that the physics 
should generically be described by two fields <j)s where s — 
a, b, with (pa = {p4't ~ l4'i)/V^ being the one entering in 
the cosine term In general, it is hard to have a complete 
form for the transformation between the (ps and the (p^. Such 
a transformation is important both for the RG analysis and 
the calculation of physical correlators which are naturally ex- 
pressed in terms of the (pa, da fields. Below, we discuss a 
special case where the transformation can be performed and 
its range of validity. 

The simplest transformation, and yet rather general, one can 
work with is a linear combination of the fields with coeffi- 
cients that are independent of the position : 

0t = Pat'/'a + Pbt'/'b ' ^ ia^9a + ib^Ob , (12) 

0; = Pai'/'a + PH'/'b ' 0^ ^ ial9a + ibiOb ■ (13) 

When p ^ q, excitations corresponding to the eigenmodes 
(pa.b carry both spin and charge modes which are respectively 
the sum and the difference of the t and | modes. As the trans- 
formation must preserve the commutation relations 

[ea{x), V<Pa' [x')] = inSaa'Six - x') , (14) 

[Mx),'^Oa'{x')] ^ i7r6aa'S{x - x') , (15) 



(11) we get that patrtaff' +pbaiba' = Saa' ■ Then, wc havc 



tat = Pbi/D , 
Ibt = -Pai/D , 



tal = -pbt/D, (16) 
ibi=Pat/D, (17) 



with the determinant 

D = patPfci - PaiPfct = (Wtfei - taihtr^ ■ (18) 

In a shortened version, we have (pa — Vcps, where V is the 
matrix of the p, and Oa = {V^^YOs- If V is unitary, the 9 and 
the (p undergo the same transformation. 

We now impose that (pa = {p(pf — q(pi)/^/2 which gives 



W^-^. (19) 



As we want to cancel the cross-terms in Eq. (fTTl . we require 
that : 



v-^K^ia^ib-^ = -v^Ki^iaiibi , 



(20) 



^PatPbt + ^PaiPbi = -9{PaiPb\ + PafPbi) , (21) 



which can be rewritten as : 

pv^fK^ibx 

Pit- +5'7)tbt 



K 



i 



qviK^ tbi = (22) 



gphbi^O. (23) 
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There exists a non-zero solution only if the condition : 



(24) 



is satisfied. When this condition is satisfied, we have a one- 
parameter family of transformations with the desirable prop- 
erty of having only one eigenmode in the argument of the co- 
sine operator The parameter is just the choice of scale of the 
field 0ft : in ID, we can change the scale of the Bose field 
provided we change accordingly its Luttinger parameter Ki,. 
Here, we choose the scale of (pb so that : 




(25) 



The condition ( l24b strongly reduces the range of applica- 
bility of the transformation : for a given coupling g, the Lut- 
tinger parameters and velocities of each species must satisfy 
the above relation. When the transformation can be used, dTTT l 
splits into a free boson field for b and a sine-Gordon model for 
a : H ^ Hoiipb) + V.sci't'a) with 'HsG = "Ho + ^2- In this 
case, the new velocities and Luttinger parameters associated 
with the a, 6 modes are given by the following relations : 



2JC 



/ 2 2 ^4. \ 

[p Wfift ^ 9 '"iKi, gpqK^Ki j , 



(26) 



velocities and Luttinger parameters then take a simple form 
in case (i), we have Va = vi, ~ vq and : 



(29) 



while in case (ii), we have K, = p^Kf and 



1 pK-t- 



9P 
2qvo' 



p^K^ 



9P 
2qvo 
3 

t 



Kf], 



In both cases, having Ka < 1 would require a very small 
(assuming q — 1 for example). This could be realized 
with long-range intra-species interactions but may not be eas- 
ily achievable. Notice a peculiarity of the formula for the mas- 
sive mode (pb in ( l29l l: while Kf and Ki are length-scale de- 
pendent (in the RG sense), the expression in ( |29] l holds on all 
length-scales. 



'Jb = 2IC ^ ^P'^'"^^^ ^ q'^v^Ki + gpqK^K{) , 



Kb = K^Ki,IC 



Vb 



where we have defined : 



/C = 



(27) 



With our definition of ipa and provided the sine-Gordon de- 
scription is applicable, the requirement for the cosine to be 
relevant, and thus to enter the single-mode phase, is simply 



Ka < 1 



(28) 



One qualitatively observes that a velocity much smaller than 
the other favors a small Ka and that large attractive interac- 
tions g < will help increase Va and reduce Ka- In the fol- 
lowing, we consider limiting cases in which the discussion 
simplifies in order to identify how the parameters would favor 
the formation of a gap in the a sector 

The limit ~ (attained with attractive interactions) sig- 
nals the transition to the phase-separated or Falicov-Kimball 
regime from the multimer phase. 



2. The limit of large asymmetry 

In order to identify the influence of the velocities ratio 
on Ka, one can introduce the dimensionless quantities v = 
vi/v^, p = q/p and 7 = g/wf - Then, (l24l) and ( l26l l are rewrit- 
ten as 



1 + p-iK^ = v{v + ^p ^K)^) , 

Ka^ + pvK\^ 



Ka = V 



— p^K^K^v + p^K^v^ 



(30) 
(31) 



If one takes into account ( [3TT i only, Ka vanishes in the limit 
of large velocity ratio — >^ or 00 and passes through a max- 
imum in between so that there are two windows of v such that 
Ka < 1. The smaller the maximum, the wider these windows 
are so, clearly, negative and large interactions (7 < 0) favor 
the mode-locking mechanism. Yet, ( l30t imposes another con- 
straint and we just consider the — > limit for simplicity. 
There, this limit is possible provided K^ ~ ~l/jp, i.e. in 
the case of attractive interaction only. As a consequence, this 
analysis shows that we should expect the formation of multi- 
mer in the attractive and large interaction regime, favored by 
large asymmetry. 



1. The case of equal velocities 



3. In the single-mode Luttinger liquid phase 



When ^ vi ^ vq, the condition ( l24l i imposes that either 



(i) g = or (ii) 



K^ 

'k[ 



The transformation and new 



Deep in the massive-a phase, one can make a crude 
quadratic approximation to the cosine operator in dS) by re- 
placing it with a mass term oc (j30^ — qcpi)"^. This leads to 
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approximate expressions for the velocity and Luttinger param- 
eter of the remaining mode b : 



2 _ P K^Vj + q K_iv-^ 
- '''''' p^K,v, + q^K,v, 



which reduces to the correct result for equal velocities. 

C. Correlation functions and the nature of the phases 

In one-dimensional models, the classification of the 
groundstates is determined by their dominant correlations. 
One can break discrete symmetries (for instance translational 
symmetry on a lattice model) but order parameters associated 
with continuous symmetries are always zero. The naming of 
a phase then corresponds to the connected equal-time corre- 
lator with the slowest decay in space. Quite generally, these 
correlators are asymptotically decaying either algebraically or 
exponentially. Such algebraic correlations are usually referred 
to as a quasi long-range order (QLRO). The slowest decay (or 
smallest exponent of algebraic correlations) criteria is based 
on a RPA argument by considering a set of weakly coupled 
Luttinger liquids [4] which shows that order will build up pro- 
vided the exponent of the correlator is smaller than two, and 
that the main instability is associated with the smallest expo- 
nent. However, if the Green's function, associated with 
which is not an order parameter, has the slowest decaying ex- 



ponent, a RG analysis shows that coupling the Luttinger liq- 
uids yield a Fermi liquid phase (provided that the decay ex- 
ponent is smaller than two again). If all physical correlators 
are exponentially decaying (apart from the density one which 
always keep, at least, a quadratic decay), the term liquid is of- 
ten used. This approach yet remains phenomenological as the 
higher-dimension situation and is much more involved. 

In this section, we follow the standard practice and con- 
sider the correlation functions of various observables to dis- 
cuss the nature of the phases that are realized in the single- 
mode and two-mode regimes. The asymptotic decay of the 
connected correlation functions associated with the order pa- 
rameter 0{x) typically reads {0{Q)0'^ {x))c oc with 
some exponent ao. In order to compute the correlators, we 
only keep the first harmonics in Eq. ^ and begin with the 
richer case of fermions where we use the representation in 
terms of right and left movers: 



(32) 



We use the results that when a field (t)a is pinned, {f{(f>a)) = 
f{{(t>a)) and its dual 6a is disordered, leading to an expo- 
nential decay. In the case of algebraic correlations, the de- 
cay exponents are obtained using the result that, for a field 
(j) described by 'Hq, the equal-time correlator associated with 
A^,n{x) = e*[™<^(^)+»«(^)l behaves asymptotically as 

(A™,„(a;)yl_„,,_„(0)) oc :r-("^^+"Vi^)/2 . (33) 
We now give the leading contributions of the order parameters 
as a function of the a, b fields, assuming general transforma- 
tion coefficients of the ^ and I modes : 



*^(.t) - e*'=-'^e-*[''<"'*''+P'"'*''-'"''^"-*''-^''l Green's function (34) 

ncr(a;) -pQcrV0a - P6(tV06 + A"^ cos(2fccrX - 2(paCT(/'a + pfca^b)) density (35) 

4't(a;)*4(a;) - e*(''^-'=^)='e"*[(P''^^P"^)'^"+(P'''^-P''t)'^''+('''-i+'''t)ea+(tH+ibt)e6l singlet pairing (36) 

- e2*t''"''^''+*'"'''''l triplet pairing (37) 



where A is a short-range cutoff. Among the multiple com- 
binations of right and left movers, we have chosen the ones 
which should lead to the lowest decay exponents, by having 
the lowest m and n constant. They usually correspond to the 
smallest wave-vector 



1. The two-mode Luttinger liquid (2M-LL) phase 

In the two mode-regime, all correlators are algebraic and 
the leading one will strongly depend on the actual coeffi- 
cients of the transformation. The expression of Eqs. (l34t- 
( [37] i are here understood with general transformation coeffi- 
cients of the t and I modes as one does not necessarily have 
to impose the restriction (fT9] l since (j)a does not here identify 



with (|9|l. The transformation coefficients can be computed 
exactly ifTsi I21I1 in the absence of the cosine term dHJ. The 
correlation functions can as well be computed directly using 
a Green's function approach ll24ll . In the presence of (|8), the 
coefficients will be renormalized in this two-mode phase to 
unknown values. This regime is rather generic and, depending 
on the interaction and densities, with many competing orders 
among which are a Fermi liquid-like phase, a superconduct- 
ing singlet or triplet FFLO phase [38] (pairing correlations 
displaying the typical — k^), a spin-density wave (SDW) 
or charge-density wave (CDW) phase. The case of equal den- 
sities, p = q = 1 has the dominant channels |4] among the 
superconducting, CDW, and SDW fluctuations. In the cases 
where spin and charge degrees of freedom separate, CDW and 
SDW states are mutually exclusive. Furthermore, for SU(2)- 
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symmetric models, x-, y- and z-components of the SDW order 
parameter are degenerate. These last remarks are no longer 
valid in our situation. 

2. The single-mode Luttinger ( IM-LL) multimer phase 

Another regime corresponds to the case where the cosine 
in Eq. ^ is relevant in the RG sense. Then, the system has 
a massive mode (j>a given by and a massless mode (pi,. 
The massless mode is described in the low-energy limit by 
a free bosonic with a velocity vi, and a Luttinger parameter 
Kf,. In this single-mode Luttinger liquid, algebraic decays 
will be ruled by this Ki, Luttinger parameter when they occur 



When the parameters of the problem satisfy Eq. (l24l . then 
the massless mode can be found explicitly. In this section 
we use these results to discuss in details the behavior of the 
correlation functions. 

When (pa gets pinned, we see that the above correlators 
(|34]|-(|37]| are all exponential because the presence of 9a in 
their expression, with the exception of the density one. In par- 
ticular, all two-body pairing channels are suppressed, even in 
the presence of attractive interactions. In order to construct 
an operator which has algebraic correlations, the prefactor in 
front of da must vanish. This is realized by taking the {p + q)- 
mer combination 'ii'^{x)'9^{x) (bound states of p 4,-fermions 
with q t-fermions) which has the prefactor qta^ + ptai which 
is clearly zero from (|T9) : 



and, in the special case of trimers, 

\i/-|^(x)\p4,(a;)*4.(a;) gik^^(,i[itat+2t^i)ea+{tbr+2tbi)0b~Pa^4'a-Pb^4'b] ^ (-39) 



in which Qqp — Sqk^ — Spki and Sp = p,p — 2, . . . , (0 or 1), 
— (?, <z — 2, . . . , (0 or 1) are integers accounting for the 
combination of left and right movers. We have used a some- 
what symbolic notation: by 'i'P{x), we mean 'i'{x + 6i)'^i^x + 
62) ■ ■ ■ '^{x + 6p), where \Si\ < A where A is the short-range 
cutoff. We stress that the family of operators d3Ft is different 
from the "polaronic" operators introduced in Ref. 18: the lat- 
ter are constructed specifically for minimizing the decay ex- 
ponents in the massless phase of ( fTTT ). On the contrary, the 
family (|38] | arises naturally in the massive phase of Eq. (fTTl l. 
as a many-body consequence of a existence of {p + (7)-body 
bound states in the microscopic counteipart of (ITll . 
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The exponent of the propagator of the {p + q)-meT then reads 
1 



^ K. 



eff 



i^eff) with the effective Luttinger parameter 



Kb 
2/C2 



(40) 



The effective theory of this (p + g)-mer object is 
then governed by the gapless mode b. Remarkably, as 
qibf + pibf = q/pbf = \/2/C, the exponent is parametrized 
only by Kb, JC and p, q. In order to have the smallest 
exponent, we have to select the combination (sg,Sp) which 
minimizes the coefficient in front of (pb (one cannot have the 
combination ppbf — qpbf = 0) and which is proportional to 
Cqp = SpP — Sqq. We list below the coefficients and cor- 
responding wave-vectors for the simplest commensurabilities: 



In this phase, the connected density correlations Af^a' (x) = 
{n„{0)n^i{x)) — (rig. (0)) (rig./ (x)) remain algebraic with the 
following dominant contributions : 



Uff{x) = 



K^^ 

' 27r2 x^ 



A 



A 



cos(2fc-|-a;) 

^2q^K,„ 

cos(2fc^a;) 



(41) 
(42) 



where Aa-a' are non-universal amplitudes. The main re- 
marks are that (i) the ratio of the zero-momentum fluctua- 
tions is exactly [q/pY while the ratio of the density is q/p 
and (ii) the wave-vectors are different since kf = 7rn— 2— and 

^ ' I p+q 
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p+q 



as well as their exponents which ratio should be 
{q/pY exactly. Notice that for the sine-Gordon model, the ra- 
tio of the amplitudes j4^-|-/A^^ are exponentially small inp—q 

Q. 

When Cqp — 1, we see that the multimer is effectively be- 
having as a spinless fermion (as expected from the combina- 
tion of a total odd number of fermions) which Fermi level is 
and Luttinger exponent Xeff- For instance, trimers belong 
to this ensemble. The effective interaction between these spin- 
less fermions, which are spatially extended objects, is highly 
non-trivial and certainly depends on the distance, density and 
microscopic parameters (a discussion of such interactions in 
the case of a boson mixture can be found in Ref. ^Ad). How- 
ever, its overall effect can be captured by K^a with effective 
repulsion expected when Xgff < 1 (dominant CDW fluctua- 
tions), and effective attraction expected if K^a > 1 (dominant 
trimer-pairing fluctuations). The latter turns out to be a super- 
fluid phase of trimers. By associating an even total number of 
fermions, one should effectively expect build a bosonic-like 
multimer Yet, we see that, in the propagator of the multi- 
mer, one cannot suppress the contribution from the 0h field 
(as Cqp ^ 0) and the exponent is not simply l/2A'eff and thus 
not simply related to the one of the density correlations as 
one would get for a simple bosonic propagator Furthermore, 
while the momentum distribution of a boson would usually 
have a peak at zero-momentum, we see that this observable 
will be here diverging at Qqp ^ 0. 



3. The case of a bosonic mixture in the single-mode phase 

As previously mentioned, the effective theory under study 
can be as well applied to the situation where the particles are 
bosons. In the single-mode phase, a bosonic multimer phase 
will emerge under the mode-coupling mechanism and the mo- 
tivation of this small section is to discuss the form of the cor- 
responding correlators. We assume repulsive interactions for 
the intra-species channels (for stability reasons and also to 
lower the K^r to be able to fulfill the Ka < 1 requirement) 
but attractive interactions in the inter-species channel (as for 
the fermions). The boson creator operators are bosonized as 
b„ ~ e*^" (dropping the higher harmonics term of Eq. (|2)) 
which immediately yields 



A(qiat+piai)Oa+iqibt+pibi)et] 



The {p + q)-meT is then a true bosonic molecule with an effec- 
tive Luttinger parameter which is exactly given by (l4Qt . The 
density correlations do not depend on the statistics and still 
have the form of (l4ni-(l42]i. 



D. Lattice commensurability effects 

So far, we have only considered two-component fluids in 
the continuum limit which is expected at generic densities on 
a lattice or in continuum space. In this section, we briefly 
discuss the additional effects arising from the presence of a 



lattice [? ]. An underlying lattice with period oq can be 
viewed as a periodic external potential, in which particles have 
a momentum being only defined modulo the reciprocal lattice 
vector 27r/ao. Therefore, umklapp processes with momentum 
transfer of a multiple of 27r/ao are allowed at low energy. If 
a Fermi momentum k^r of a species cr, is itself a multiple of 
27r/ao, i.e. if a density of species a is commensurate with 
the lattice, sria — integer/ao, with an integer s, an additional 
term cos(2s0ct) appears in the low-energy Hamiltonian. The 
effects stemming from such a cosine operator alone are well 
known : for K^r ^ 2/s^ the cosine is relevant in the RG sense 
and the system undergoes a Mott transition into a density wave 
state with the unit cell of s lattice sites. In a two-component 
system, it is possible to have two operators of this sort, one 
for each species. Furthermore, if the densities are such that 
sn^ + s'ni is an integer (we set ag = 1 from now on) for 
some integers s and s', there is yet another term in the low- 
energy Hamiltonian, namely cos 2(s(/)^ + s'4>i) (cf. Eq. 
Here, we analyze a simple special case where 



p - g nj, = , 
r 71^ + ; nj, = 1 , 



(43) 
(44) 



or n-j- = q/ {pi + qr) and — p/{pl + qr) with the integers 
p, q, r and I. Given (|43]l-(l44ll. Eq. ([3]) yields the Hamiltonian 
in the form Ho(?!'t) + ^o(04.) + Hi + 'Hcos with 



He 



Gi / dx cos2(p0^ -q(t)i) 



+ G2 J dx cos 2(r + I 04_) 
+ G3 dx cos 2{pl + qr) (f>^ 



+ G4 dx cos 2{pl + qr) (j)^ 



(45) 



(46) 



(47) 



(48) 



where Gi,. .. ,Gi are non-universal amplitudes. Interpretation 
of Eqs. (|45]|-(|48]| is straightforward: Eq. (l45t stems from the 
condition ( |43]) and is thus insentive to the presence of the lat- 
tice (cf. Sec. Ill Al ): Eqs. (gTll and (|48]) favor the Mott localiza- 
tion of the species \. and respectively. On the other hand, 
operator (l46t is unique to two-component lattice systems and 
owes its existence to the peculiar commensurability condition 
(|44| ). The physical meaning of ( l44l l is clear: by analogy with 
Sec. Ill CI it favors the quasi long-range ordering of the o^sr&- 

tora+i-*^!*!)'. 

In the following, for sake of simplicity, we assume equal 
velocities of the two-components and drop the "Hi term. The 
dominant instability of the massless theory 'Hol^t) + Ho(04,) 
is due to the operator with largest positive scaling dimension. 
Depending on the values of and K^, the following inequal- 
ities define which of the operators (I45t-(l48t is relevant: 



p'K^ + q^K^ 2 , 
r^K^ + PK^ ^ 2 , 
[pi + qr)^K^ sC 2 , 
{pi + qr)^K^ 2 , 



(49) 
(50) 
(51) 
(52) 
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FIG. 2. (Color online). Diagram showing the effect of commensurate 
densities (see text for discussion) in the special case of = 1/3, 
= 2/3, i.e. p = 2 and g = r = ? = 1. For 2 - 2/9 < 
Ki < 2/9 the interaction with the lattice leads to a formation of a 
'trimer crystal' state. For larger(smaller) values of the system 
has a phase transition from a massless phase into a Mott insulator of 
t(i)-component. The 'trimer' operator cos 2(20^ — (pi) is always 
subdominant. 



respectively. In Figs. |2]and|3]we plot the {K^,Ki) diagrams 
corresponding to Eqs. (|49]l-(l52]i for two values of the den- 
sities. We see that which instability takes place depends on 
the values of the bare Luttinger parameters and K^, and 
thus on microscopic details of an underlying lattice model. 
Numerically, a crystal phase has been reported ||4T1| in a two- 
component bosonic Hubbard model and a similar result is pre- 
sented in the fermionic counterpart in Sec. IIIIFI for the com- 
mensurabilities discussed in Fig. |2] These phases do corre- 
spond to the locking of several combinations of the modes ac- 
cording to Eqs. (l45]l-(l48Tl but they are achieved for very large 
asymmetry. Consequently, the above criteria (|49^-(|52t deter- 
mined for equal velocities are not directly applicable in these 
situations. The quantitative predictions of (|49]l-(l52]i could be 
relevant to the case of strongly renormalized K^, for instance 
with long-range intra-species interactions. 

A striking feature of the phase diagrams |2]and|3]is the ap- 
pearance of the multicritical points where several instabili- 
ties compete. In the above treatment we have only consid- 
ered an effect of various operators (I45t -(|48^ alone. An inter- 
play between different operators is non-trivial and may lead 
to consequences not captured by the simple power counting 
of Eqs. (|49]l-(l52]i. Hence, applicability of the above analy- 
sis in the vicinities of the multicritical points is not granted. 
There are several possible scenarios of the phase transitions 
at such multicritical points. For one thing, it is easy to con- 
struct fine-tuned theories where two continuous transitions oc- 
cur simultaneously. An other possibility is a first order tran- 
sition, as been observed in numerical simulations of higher- 
dimensional bosonic systems |42]. Detailed analysis of these 
multicritical points is beyond the scope of the present paper 
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FIG. 3. (Color online). Same as Fig.|2]for = 1/5, nj, = 2/5. In 
this case, Eqs. I l43t -l l44t allow two sets of solutions: (A) I — 1 and 
r = 3, and (B) / — 2 and r = 1, with p — 2 and g = 1 in both 
cases. Solution (A) is always subdominant, while (B) dominates in 
the window 2/5 < < 2/25. For 48/25 < < 2/5, the 
dominant instability is the formation of a Luttinger liquid of trimers. 



III. TRIMER FORMATION IN THE ID ASYMMETRIC 
HUBBARD MODEL 

In this second part, we study the emergence of a trimer 
phase on a particular microscopic model : the ID asymmetric 
attractive Hubbard model. After defining the model and pro- 
viding its phase diagram as a function of the parameters, we 
discuss some Umitations of the bosonization approach to this 
model and an alternative phenomenological description that 
completes the interpretation of the obtained data. 

A. Model and qualitative aspects 

We consider two species of fermions which internal degree 
of freedom is denoted by a spin index a. They hop on a lat- 
tice with a spin-dependent amplitudes tc (which would ex- 
perimentally corresponds to different optical lattices for each 
species) and interact locally only in the inter-species channel 
with a Hubbard term U which we take negative, as suggested 
by the arguments of Sec. Ulland as a natural choice to favor 
bonding between particles. The Hamiltonian is then : 

■H = - ^ t„[cl_^^^^Ci^a + h.c.] + U^ni^^ni^i . (53) 

One of the key parameter for the physics is the ratio be- 
tween the hoppings rj — ti/tf. In order to have the possibil- 
ity of forming trimers, we take the commensurate condition 
2n^ = ni but the total density n varies freely and is another 
important parameter of the physics. Using the notations of 
Secini we thus have p ^ 2 and q = 1 (the simplest new com- 
bination one can have). The Fermi momenta are ka- = Tm„ 
and free fermions Fermi velocities read v„ = 2fcr sin(7rncr). 
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FIG. 4. Opening of the trimer gap increasing mass asymmetry (low- 
ering rj — ti/t^) for a fixed interaction and density. Tiie magnitude 
of the gap (in units of t^) is small in comparison to U and t^. The 
grey areas are estimates of the transition points. Inset: finite size ex- 
trapolations of the gap. The upper dashed curve shows the behavior 
for r; = 0.2 when entering in the FK regime. 




FIG. 5. (Color online). Examples of fits of the entanglement entropy 
in the two-mode and single-mode phases using Eq. l l55t . It shows a 
clear quantitative difference with respectively c ~ 2 and c ~ 1, as 
expected. 



B. Phase diagrams 



Since nj^ < 1, the maximum total density one can have for 
this commensurability is n = 3/2. 

The above Hamiltonian has been widely studied in the case 
of balanced [3] and imbalanced densities |j43J but the spe- 
cial commensurability where trimers emerge has only been 
investigated for one set of data in Ref. showing that the 
pairing correlations were indeed suppressed, in agreement 
with the bosonization approach. When the asymmetry is very 
large, one species behaves quasi-classically (they get local- 
ized) and the model is in the regime of the Falicov-Kimball 
(FK) model | 44] where there exists a lot of quasi-degenerate 
states at low-energies, analogue to a phase separation regime. 
We expect generically a first-order transition to this segregated 
(or demixed) phase when lowering 77 in the phase diagrams. 
The FK regime can display rather rich physics recently inves- 
tigated in Ref. |45| and which will not be analyzed here : our 
aim is only to draw the boundary of this regime. Numeri- 
cally, the transition to the FK is rather sharp and all observ- 
ables clearly display segregation. For -q — 1, the arguments of 
Sec. HI] suggest that the two-mode regime will be generically 
realized. Qualitatively, in a strong-coupling picture where two 
spin-| fermions are localized on neighboring sites, the delo- 
calization of a spin-^ electron on these sites will be favored by 
attractive interactions, forming a very local trimer state. This 
picture will be correct at small enough densities and actually 
not too large U and too small ry, otherwise such bound states 
will agglomerate with other spin-f and \. fermions, leading 
to the FK regime. We thus expect the formation of the trimer 
phase in the vicinity of the FK but at both finite U and finite 77. 
Within the framework of Sec. |II]and considering that the start- 
ing point of bosonization are free fermions, the ratio between 
the velocities v ~ V]^/v^ — 2rjC0&{'Kn/'i) supports that small 
■q clearly favors the formation of trimers while small densities 
should not. 



The phase diagrams of model ( l53T l are numerically deter- 
mined using standard DMRG with open-boundary conditions 
(OBC) and keeping up to M = 2000 states. In order to dis- 
criminate between the different possible regimes, we use both 
"global" probes and local observables and correlation func- 
tions. Among "global" probes, one can use the trimer gap At 
associated with the formation of the bound state. It can be 
defined following Ref. 31 as: 

A, =Eo{N^ + 1, iV^ + 1) + Eo{N^, + 1) 
-Eo{N^ + l,N_^ + 2)-EoiN^,Ni) , 

with Eq{N^,N^) the ground-state energy with N^,N^ 
fermions. Results as the function of the asymmetry 77 for 
an incommensurate density n = 3/7 and large interaction 
U = — 4t^ have been extrapolated to the thermodynamical 
limit and are given in Fig.|4] The slow opening of the trimer 
gap is qualitatively compatible with the sine-Gordon behavior 
of Sec. although the transformation is not directly applica- 
ble for any 77. Notice that the whole system remains gapless. 
The slow opening of the gap makes it difficult to precisely lo- 
cate the transition point. In such situation, a usual approach 
would be to use the prediction on the critical Luttinger pa- 
rameter Ka = Ka '^he transition point. Furthermore, the 
determination of Ka using coiTelators in the two-mode phase 
is very difficult as it would require to know, and then to dis- 
entangle, the complicated expression of the exponents as a 
function of Ka and Kb to extract them independently. 

Therefore, we use another global approach to the distinc- 
tion between the two-mode and single-mode phases, which is 
particularly well-suited for this model, and more generally in 
a similar context. Using universal results on the entanglement 
entropy (EE), the central charge c of the model can be ex- 
tracted which directly gives access to the number of bosonic 
modes, without further information on their nature. Hence, 
we expect c = 2 in the two-mode regime while c = 1 in the 
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FIG. 6. (Color online). Central charge c obtained from fits as in Fig.|5] 
as a function of the asymmetry for 17 = —4tf and n — 3/7. The 
stair-like behavior with increasing system size L allows an efficient 
determination of the transition. In the Falicov-Kimball regime (left), 
fits give c ~ or irrelevant numbers, (a) Zoom of the c(L) curves 
in the crossing region illustrating the extraction of crossing points 
between successive sizes, (b) Tentative finite-size extrapolation of 
the crossing points Tjc{L). 

single-mode trimer phase. This stair-like expectation in the 
thermodynamical limit will be smoothed out by finite-size ef- 
fects. The central charge is obtained on finite-systems using 
the following ansatz for the EE between a left block of size a; 
and the right block of length L + 1 ~ x with OBC: 

S (x)^- In d( x\L + 1)+ At( x) + B , (55) 
6 

where d{x\L) is the cord function 

d(x\L) = - sin (—) , (56) 

TT \ L / 

and t{x) is the local kinetic energy on bound (x, x + 1) (ob- 
tained numerically), and A, B are fitting parameters. The first 
log term is the leading and universal one |46] while the sec- 
ond accounts for finite-size oscillations due to OBC and which 
can have a significant magnitude ll34ll47ll . It is thus essential 
to take them into account to improve the quality of the fits. 
In the end, there are only three parameters in the procedure 
and typical examples in both the two-mode and single-mode 
phases are given in Fig. |5] Systematic fits on finite-size sys- 
tems provide an estimate of c as a function of the parame- 
ters. As seen in Fig. |6] the c{L) curves crosses around the 
transition point. Although we do not have any quantitative 
prediction for the finite-size corrections of c(L) obtained in 
this way, we can argue that if L is smaller than the coiTela- 
tion length associated with the tiimer gap, c{L) will be larger 
than one as the system is effectively in a two-mode regime. 
Thus, c{L) should decrease with L towards one in the single- 
mode phase, as observed. In the two-mode phase, there is no 
obvious discussion : we only expect that the larger the sys- 
tem, the better the agreement with the continuous limit. One 
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FIG. 7. (Color online). Maps of the central charge c vs interaction 
U and asymmetry r] for a system with L = 112 at four different 
densities. For n = 3/7 ~ 0.428, the lines with error bars are the 
ones estimated from Figs. |4] and |6] The rj = Q cuts correspond to 
data obtained with a very low but non-zero value r; = 0.005. 



can also check the effect of the number of kept states on the 
fits and see that they don't have the dominant effect in this 
model which converges well numerically. We have estimated 
the transition point by extrapolating the crossing points be- 
tween successive sizes (see Fig.|6ta)) as a function of the in- 
verse size (see Fig.|6lb)). From this approach and the opening 
of the gap, we get a critical value rjc ~ 0.54 ± 0.02 for the 
mass asymmetry on this cut. Although the gap is rather small, 
the trimer region appears to be rather wide. 

Using the central charge calculation, one can map out the 
phase diagram in the (77, U) plane for a fixed density, or in the 
(77, n) plane for a fixed interaction U . Results are gathered in 
Fig.|2]and[8]respectively. These diagrams display bare data for 
a given system with a rather large size L = 1\2 and the previ- 
ous estimate of the cut is given as eiTor bars. These diagrams 
show that a wide trimer phase can be achieved at large enough 
interactions, small enough 77, as expected, and also that low 
densities strongly favors their formation. At large densities 
71 ~ 1.3, the tiimer region vanishes within our grid resolution 
so that it is at most confined to a very tiny region between the 
two-mode phase and the FK regime. While the large-|[/| situ- 
ation is rather clear, the competition between the three regimes 
at small U is more involved. Indeed, two scenaiios can occur 
in the (7/, U) plane: either the trimer phase always separates 
the FK and two-mode regimes, corresponding to two bound- 
aries starting from the (77 = 0, C/ = 0) corner, or there is 
a critical \U\ above which the trimer phase emerges, corre- 
sponding to a triciitical point [rjc, Uc)- We could not numer- 
ically discriminate between both scenaiios, but we do find a 
small trimer region at relatively small [/s ([/ ~ —1 —2 1^) 
for most densities : we do not have evidence for a tricriti- 
cal point with a large Uc- As the density plays a central role 
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FIG. 8. Map of the central charge c vs asymmetry ri and total density 
n for fixed interaction U = — 4f-f on a system with L — 128. The 
lines with error bars are the ones estimated from Figs.2]and|6] 

in the stabilization of the trimer phase, we give in Fig. |8]the 
central charge map for a fixed interaction U = —At-^ as a 
function of the total density n and mass asymmetry. A similar 
question about a intervening trimer phase between the two- 
mode and the FK regimes can be raised. While the two-mode 
and FK are clearly separated at small densities, we found that 
if a trimer intermediate phase exists at large densities up to 
the {rj — 0,n = 1.5) point, its extension will be particu- 
larly small (not seen within our numerical calculations). In 
addition to the three main phases, commensurability effects 
are also present in this diagram. When the maximum den- 
sity n = 1.5 is reached, the |-band is completely filled while 
the f-band is half-filled, leading to a single-mode phase well- 
captured by the central charge approach. Lastly, as it will be 
discussed in Sec. 1111 Fl a crystal phase (fully gapped) exists 
for the commensurate density n = 1 at very small rj and is 
pointed on Fig. [8] Other commensurabilities could yield ad- 
ditional crystal-like phases in this diagram but this is beyond 
the scope of this study. 



C. Observables and effective behavior of the trimer liquid 

In this section, we give the behavior of several observables 
in order to see how they are affected by the entrance into the 
trimer phase or the FK regime and, as well, to investigate the 
effective behavior of the trimer fermion. 



1. Local observables 

First, we select a set of local correlators (living on sites or 
on bonds) which illustrate the phenomenological picture of 
the different parts of the phase diagram. We compute the lo- 
cal double occupancy {n,i^^ni,i), the trimer local operator as 
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FIG. 9. (Color online). Maps of averaged local quantities in a system 
with L = 56 at total density n — 3/7 ~ 0.428. The non-interacting 
expectation ((7 = line) has been subtracted in order to unveil the 
effect of the interaction. 

7i = i((ni,t«i, 4, '^i+i.;) + ("-»+i,t'ii,;"-i+i,i)) (since the light 
particle is in principle delocalized above two heavier), and the 
density correlators {ni^^ni^i,-^) and (n-i^nj+i^^). This choice 
of local correlators is well suited to a strong-coupling picture 
as pairs or trimers should in principle correspond to a narrow 
bound-state, spread over only a few lattice sites. These local 
correlators should then pick up a reasonable weight of the lo- 
cal bound-state. The results are averaged over all lattice sites 
and plotted in Fig. |9l The expectation value at U — has 
been subtracted so that the reference state is the free fermions 
limit at a given 77 (the pairing or trimer local correlators de- 
fined above are obviously non-zero even in the free fermions 
limit). Fig.|9]display behaviors in qualitative agreement with 
the picture we have on the trimer formation : the density 
correlator is nearly zero everywhere but in the FK regime, sig- 
naling phase separation. On the contrary H density correla- 
tor increases significantly in the region corresponding to the 
trimer phase, surrounding the FK pocket, and together with 
the local trimer density %■ Lastly, we see that the double oc- 
cupancies, or pairs, acquire a strong weight with negative U 
everywhere in two-mode and single-mode regions : they are 
either "independent" or embedded in the trimer bound-state. 
Their coherence can yet be probed only by measuring corre- 
lations as discussed below. 



2. Pairing and trimer correlations, effective behavior of the 
trimers 

We now turn to the behavior of correlation functions across 
the phase diagram. From Sec. Ill CI and as already observed 
for a particular point in Refs. ^ andi31j the pairing correla- 
tions change from algebraic to exponential decay when enter- 
ing into the trimer phase. These correlations are here com- 
puted in the singlet channel and for local pairs Pi = Ci^^Ci j,. 
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FIG. 10. (Color online). Behavior of (a) pairing and (b) trimer corre- 
lations when lowering rj and entering the trimer phase along a cut at 
U — —4t-f in the phase diagram (absolute values are displayed). The 
Inset of (a) shows the same data but in log-linear scale to highlight 
the exponential decay. 



In addition, we compute the trimer correlator using the local 
trimer operator Ti = Ci^^Ci^iCi^i^^ defined on neighboring 
sites. The associated correlation functions = {PjPi^x) 
and T{x) — {TjTi+x) are computed with i taken at the cen- 
ter of the chain. Increasing the mass asymmetry along the 
same cut aiU — —At-\ as in previous figures, the suppres- 
sion of pairing correlations is clearly seen in Fig. [TOl a). On 
the contrary, trimer correlations, which are subdominant in 
the two-mode regime, are boosted by smaller 77, both in am- 
plitude (as for the local correlators previously evoked) and in 
the decay exponent which gets smaller (see Fig. [TOl^b)). No- 
tice that the wave-vector is the same for both correlators since 
/cj.— /cf = /cf = 7rn/3 for this commensurability. We have ten- 
tatively extracted the decay exponents of both correlators by 
fitting the functions using a power-law modulated by a cosine 
oscillations. The correlation length ^ of the pairing correlator 
in the trimer phase is obtained using an exponential envelope 
g-a;/C -pjjg results are gathered in Fig. [TT] showing the evo- 
lution in both phases. We must stress that the data computed 
on a finite-size system display a transition at a lower rj than 
in the thermodynamical limit. From Sec. Ill CI we expect that 
the decay exponent of the tiimer propagator is of the form 
(i^eif + K^tt)/'^ while the f-density correlations have a de- 
cay exponent of 2K^{i. A first consequence is that the trimer 
exponent should always be larger than one which is not re- 
produced for the lowest 77s, and which we attribute to numer- 
ical inaccuracies of the fits of trimer correlations. Besides, 
inverting [K^fi + K~^)/2 to get Xeff is subjected to strong 
errors when K^s ~ 1 and does not tell whether Xetf > 1 or 
ifeff < 1 which is essential for the effective behavior of the 
trimers. In order to get a better estimate of Xetf, we rather 
use Friedel oscillations of the ^-density operator which decay 
exponent a is equal to A'eft in the trimer phase according to 
Sec. Ill CI Even though the approach of Sec. Ill CI is not ap- 
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FIG. 11. (Color online). Comparison of the decay exponents of pair- 
ing and trimer correlations. The correlation length of pairing correla- 
tion (in units of the lattice spacing) is given in the trimer phase. The 
strong grey areas are the previous estimates of the transition points 
while the light maroon area illustrates the location of the transition 
on the L — 147 finite-size system under study. The exponent of 
the Friedel oscillations of Uf is also displayed, together with the ex- 
pected trimer exponent derived from it (see text for discussion). 



plicable for most parameters, the fact that there exists an ef- 
fective Luttinger exponent K^ff describing the physics of the 
fermionic trimer with a propagator (ifeff + K^^)/2 and with 
Friedel oscillations with Xeff is more general : the limitation 
of the bosonization approach is rather that K^ff will not take 
the form of Eq. (|40] |. Local observables are believed to have 
less numerical errors associated with a finite number of kept 
states than correlations ll48ll . Thus, we fit the Friedel oscilla- 
tions of the t-density using the following symmetric ansatz 



cos q{i - ^) 



(57) 



with 1 < X < L and only four fitting parameters no, A, q 
and a [? ]. In the trimer phase, we expect a — K^ff. Some 
typical fits are plotted in Fig. [T2l a). From them, we extract 
the decay exponent a and plot it on Fig. fT^l b) as a function 
of the total density. A cusp is found around n ~ 0.6 signaling 
the transition from the two-mode regime to the trimer phase. 
We have seen that a low density favors the formation of the 
trimer phase. This figure shows that, in the trimer phase, we 
have Keff < 1 for the larger densities, corresponding to a re- 
pulsive effective interaction between trimers (dominant CDW 
order of trimers). We observe that the exponent increases with 
decreasing density, compatible with the fact that at low den- 
sities in the trimer regime, the trimers should be close to free 
spinless fermions having Kf^g ~ 1. Bare data for the small- 
est density on a system with L ~ 144 even display an expo- 
nent A'eff ~ 1.2 larger than one. Interestingly, the trimers in 
this model are necessarily objects with a finite extension of at 
least two sites and two close trimers may have the possibility 
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FIG. 12. (Color online), (a) Typical Friedel oscillations of the t- 
density for U — —4t-f, rj — 0.3 and L — 144 for various densities 
n. Full lines are fits using Eq. ( I57t . (b) Decay exponents obtained 
from the fits as a function of the density. The cusp at n ~ 0.6 roughly 
corresponds to the transition from the two-mode to the single-mode 
regime, (c) Large exponents at low densities, close to the FK regime 
when lowering ry : increasing the size tends to reduce the exponent 
below one. 



to overlap by delocalizing their t electrons. The distance de- 
pendence and sign of the effective interaction between trimers 
is non-trivial - a perturbation theory to derive it looks chal- 
lenging as it involves many sites and degrees of freedom. Yet, 
since the trimer phase is found close to the FK regime, we 
can expect the effective interaction to become attractive close 
to this boundary, leading to ifeff > 1- Such physics would 
coiTespond to a superfluid phase of trimers. This would be 
physically very remarkable since the microscopic Hamilto- 
nian (ISJt would contain both the formation of bound-states, or 
molecules, and their effective superfluid behavior However, 
the behavior close to a phase separation at small densities is 
numerically involved. Indeed, increasing the system size L 
shows that a actually tends to decrease below one, as reported 
in Fig. fTSl' c). or one enters the FK regime for larger sizes. We 
did not find clear evidence of a stabilization of K^-ff > 1 in the 
thermodynamical limit and understand the observed ifeff > 1 
as finite size effects. A superfluid droplet picture can be qual- 
itatively put forward. Starting from the FK regime and look- 
ing at the local density pattern, one sees that the fermions 
are clustered into droplets while other parts of the box are 
empty. Approaching the trimer phase from the FK regime 
tends to increase the size of these droplets to gain kinetic en- 
ergy. When a box confinement is present (finite system with 
OBC), it naturally favors the overlap between trimers, by de- 
pleting the edges, and can prevent the droplet to form (for 
instance if their typical size is larger than the box size). In- 
creasing further the size of the box (at constant density) can 
lead to droplets formation. This is a possible interpretation of 
the data observed in Fig. [T2l c). In addition, we must stress 
that there is a lot of competing low-energy states in the FK 
regimes so DMRG, as an essentially variational methods, can 
be trapped into metastable states. Even though the thermo- 
dynamical limit is unclear, it is experimentally motivating to 
have signatures of superfluidity on mesoscopic confined sys- 
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FIG. 13. (Color online). Density correlations in (a) the trimer phase 
and (b) the two-mode phase with n = 3/7. 



tems as one has in cold-atoms setups. We further mention 
that a recent careful study of the t-J model on a chain [491 
which could qualitatively contain a similar phenomenon as 
pair-clustering did not find evidences for such clusteiing. 



3. Locking between density correlations 

Lastly, the comparison of density-density correlations in the 
t and i channels is another interesting point of this model. In 
fact, the bosonization results of Sec. IllCl predicts that the ex- 
ponent of A/f-|- (x) should be four times larger than the expo- 
nent of A/4,4, (x) (if both remains smaller than two), and that 
the dominant wave-vector should differ by a factor two. Nu- 
merically, the typical behavior for a rather large interaction 
U = —At^ is given on Fig. [13] for two values of the asym- 
metry f] in the single-mode and two-mode regimes. We see 
that in the two-mode phase (Fig. [T3l b)). the two fluctuations 
have a slightly different exponent, that the amplitude are quite 
different (including the natural factor four). Yet, the domi- 
nant wave-vectors are both 2k^ = k^. In the trimer phase, 
the disagreement with the bosonization picture is even worse 
since the two densities are locked together, and nearly iden- 
tical (Fig. [T3l a)). This latter fact cannot be explained by 
the l/x^ decay since the leading term is the oscillating one, 
with an exponent clearly smaller than two. It is yet physi- 
cally not surpiising in the strong coupling picture of Fig. [T]: 
trimers are local bound-states separated by the typical distance 
l/rif — 2TT/2k^ which does correspond to the 2k^ fluctua- 
tions but cannot be accounted by any of the harmonics for the 
^-density operator of Eq. (|3]l (we work at an incommensurate 
filling). This short-distance binding cannot be captured by the 
bosonization results of Sec. |II]but a phenomenological Bose- 
Fermi approach described in the next section can account for 
this strong-coupling regime. Lastly, the same comment can 
be made on Friedel oscillations on the ^-component : they are 
locked to the ^-componentin the strong-coupling picture. One 
might argue that there could be a crossover from the weak- 
coupling to strong-coupling picture of Fig.[T]as \U\ increases. 
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so that the bosonization results could be valid in the small Us 
region. However, the trimer region is very sharp at small Us 
and we could not find evidence for such a weak-coupling be- 
havior in our numerical data, although we cannot exclude such 
a possibihty. 



picture, the trimer is simply a bound-state between the bosons 
and the fermions so its propagator reads 



(62) 



^ikfX^i[{idf+idb)6d + iiaf+icb)6a~Pdf(f'd-Pcf<t>c] (-^2) 



D. Phenomenological Bose-Fermi picture at large \U\ 

We here propose a simple picture that reconciles the nu- 
merical observation and a bosonization approach at the prize 
of a strong assumption, physically reasonable at large neg- 
ative U, but difficult to justify rigorously starting from the 
microscopic model. This picture has been for instance pro- 
posed at large interaction and low-density limit [50]. Bose- 
Fermi mixtures of ID models have extensively studied in re- 
cent years |5 ll - l54ll and a similar 
picture emerges in certain regimes of three-component Fermi- 
gases |55J. When \U\ is large, t and | fermions naturally form 
onsite pairs which are effectively hard-core bosons which we 
label b. We phenomenologically assume that the system is 
equivalent to a Luttinger liquid of hard-core bosons with den- 
sity Tib — n^, an effective velocity vi, and Luttinger parameter 
Kf,, while the remaining unpaired I fermions behave as a Lut- 
tinger liquid of fermions labeled by / and with parameters 
nf,vf and Kf. These two Luttinger liquids interact through 
an effective interaction which will contain terms such as 



dx cos [27r(ny — ni,)x — 2(0^ — 0^)] 



(58) 



which have the tendency to lock the fields and (j>i, to- 
gether (with {(j>f) — {(j>i,) for attractive interaction), provided 
n-f — nt- Such an effect has already been discussed in the 
context of Bose-Fermi mixtures ||53[ l54ll . Clearly, the latter 
relation is the same as the trimer commensurability condition 
ni = 2n^ (because bosons carry two particles) so that the for- 
mation of trimer is now interpreted as a bound-state between 
the bosons and the fermions. Following the same reasoning as 
in Sec. Ill B I we can introduce a general transformation of the 
b/ f fields into two new fields c/ d where (pc — {<j) f — (pb) / \/2- 
Writing the matrix transformation from the r = c, d to the 
s = 6, / as prs for the 0s and trs for the dual 9s, we have 



1 



V2 



icb = — 



1 

72' 



(59) 



which is slightly different from Eq. ([19). Yet, the 2fci?-like 
fluctuating part of the density correlators for the fermions and 
the bosons will have the leading contributions (dropping the 



terms) : 



A/'//(^) 



cos{2kfx) 

X * 
cos(2fc;,a;) 



(60) 
(61) 



which have both the same wave-vector associated with the 
Fermi levels kf = kb = and same decay exponents since 
pdf — Pdb from the canonical transformation relations. In this 



From Eq. i59i and the determinant of the transformation ma- 
trix, which gives that tdf + tdb — — 1/P<#' '^^ obtain 
that the propagator is of the spinless fermionic type with an 
effective Luttinger parameter K'^g = p^j:Kd- Clearly, both 
the /-fermions and 6-bosons propagators become short-range, 
the latter corresponding to the pairing correlations in the na- 
tive fermionic model. Consequently, we recover the physics 
of the trimer phase developed in Sec. |II] with a better agree- 
ment with the numerical observations in the strong coupling 
regime. However, splitting the initial gas of I fermions into 
two parts can only be done phenomenologically and could be 
questionable in a microscopic derivation. This highlights the 
limitation of the bosonization approach of Sec.|II]at short dis- 
tances (high energies). 



E. Possible observation of the trimer phase in the presence of 
parabolic confinement 

In this section, we briefly discuss the condition to favor the 
trimer phase in the presence of a parabolic confinement, as 
used in cold-atoms experiments. Our goal is only to exhibit 
some parameters for which the trimer phase is stabilized and 
to give some qualitative comments. The trapping potential is 
taken into account by adding the quadratic term 



U 



1 



trap 



(64) 



to Eq. |53] with the trapping frequency uj and the center of 
the lattice io ~ {L + l)/2. According to a local-density ap- 
proximation (LDA) picture and using the phase diagram of 
Fig. [8] the trimer phase is likely to be found at small enough 
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FIG. 14. (Color online). Realization of the trimer phase in a trapped 
configuration with parameters U = —it^, i] = 0.3, iu = 0.002-^17' 
and A'^ = 51 fermions. (a) Local densities profiles, (b) Correlation 
functions from the center of the trap. Inset: same in log-linear scale. 
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FIG. 15. (Color online). Observation of the crystallization of trimers 
in the commensurability situation of Fig.|2](n = 1). (a) the short- 
range behavior of both pairing and trimer correlations, (b) ordering 
of the local density with the expected period of three sites. 



densities and not too small rj to prevent the occurrence of the 
FK regime. However, we find that the average density of the 
trapped system is strongly dependent on the Hamiltonian pa- 
rameters. At fixed number of particles N and trap size lu, 
changing U and 77 strongly affects the radius of the cloud and 
the density in the middle. We only exhibit in Fig. [T4lparame- 
ters for which the main features of the trimer phase are repro- 
duced in the presence of a parabolic confinement. The den- 
sity profiles of Fig. [T4l a) illustrate the locking of the t and | 
densities (up to exactly a factor two), and the emergence of 
an appreciable density of local trimers (each local maximum 
roughly corresponding to a trimer). In Fig. [l4T b). the pairing 
and trimer correlations are strongly different from that of a 
superfluid phase : we have dominating trimer correlations and 
exponential pairing correlations as in the homogeneous coun- 
terpart. In agreement with a LDA picture, since we have seen 
that K^s decreases with density, the trimer correlations are 
boosted at long distances. Similarly, the pairing correlations 
decrease slightly faster then an exponential close to the edge 
of the cloud. These results are encouraging in the perspective 
of a possible achievement of the trimer phase in actual exper- 
iments. 



F. Observation of a crystal of trimers 



According to the analysis of Sec. Ill Dl a crystalline phase of 
trimers can occur in this lattice model when the total density 
n is commensurate. Evidence of this scenario together with a 
phase diagram for n — 1 has been proposed in Ref. 41 in the 
case of a mixture of two-component bosons, for large enough 
asymmetry. As the order parameter (the density) associated 
with this transition is independent of the statistics, we expect 
a similar scenario (see Sec. lllDI l and a similar location of the 
transition in the fermionic version of the model under study. 
Indeed, we give in Fig. [15] an example of the crystal phase. 
Notice that very small 77 are required to stabilize such a phase. 
We have not investigated the extension of the phase which 
should rather be small on the scale of the phase diagram of 



Fig. [8] and its neighboring phases which could be either the 
two-mode LL or the trimer phase. A crude argument can be 
proposed to understand this crystallization within the Bose- 
Fermi picture of Sec. 1111 Dl : when the mass asymmetry is very 
large (very small rj), it is reasonable to say that the mass of 
the boson will be essentially the one of the heaviest particle, 
which is the same as the unpaired fermions so that ~ w;,. In 
terms of commensurability effects, one has nj = rib = n/i 
so that standard umklapp terms at 27r(n/ + rib) — 47r/3 do 
not account for the crystallization. One rather has to look for 
higher order terms with commensurabilities such as 2n/ + 
rib = Uf + 2nb = n = 1, which are typically associated with 
terms like 



dx cos2{2(f>f + (f)b) , 



(65) 



in addition to the one of Eq. ( 1581 ). Such term can lock the field 
2<j)f + 4)b and make the system fully gapped. Lastly, we would 
like to stress that such commensurabilities are rather surpris- 
ing in terms of the initial fermion densities as they belong to 
odd filling fractions — 1/3 and — 2/3. 



IV. CONCLUSIONS 

In summary, the consideration of unusual commensura- 
bility conditions in density-density interactions for ID two- 
component gases leads to a very rich physics with the possi- 
bility of building bound-states of [p + -particles as the lead- 
ing order Such a mode-locking mechanism can be desciibed 
within the framework of Luttinger liquid theory which reveals 
the main ingredients to stabilize such a new phase. In partic- 
ular, mass or velocity asymmetry is shown to diive efficiently 
the transition into the multimer phase. Novel fully gapped 
phases are proposed when taking into account umklapp cou- 
plings specific to lattice models at commensurate densities. 
These ideas are illustrated and confronted with the asymmet- 
ric ID attractive Hubbard model for the special commensura- 
bility 2n^ = for which the formation of trimers is found. 
The features of the phase diagram are computed, display- 
ing the important role of the density in favoring the trimer 
phase. The effective behavior of trimers, which are effec- 
tively spinless fermionic objects, is very sensitive to the den- 
sity and mass asymmetry. Although the model seems to have 
promising features to sustain a superfluid phase of tiimers, 
we did not find clear evidence for it in the thermodynamical 
limit, while finite-size systems display a "superfluid droplets" 
physics. Notice that superfluidity of bound-states made of 
four fermions (quartets) can be achieved reliably in ID with 
a four-color Hubbard model ll33ll34il . There, the bound-states 
are bosons which natural "free" regime (attained in the low- 
density limit) is a superfluid phase. A superfluid phase of 
trimers would in this respect be even more exotic but is in 
strong competition with phase separation. Lastly, we found 
that a trapping confinement supports the trimer phase for rea- 
sonably high densities and that surprising crystal phases can 
emerge at commensurate densities. 
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